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We analyse a nonadiabatic self-consistent field method by means of an exactly-solvable model. The 
method is based on nuclear and electronic orbitals that are functions of the cartesian coordinates in 



o 

q \ the laboratory-fixed frame. The kinetic energy of the center of mass is subtracted from the molecular 

Hamiltonian operator in the variational process. The results for the simple model are remarkably 
accurate and show that the integration over the redundant cartesian coordinates leads to couplings 
among the internal ones. Present results show that the nonadiabatic SCF method in cartesian 
coordinates does not always yield upper bounds because the variational principle is violated. 
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^ ' I. INTRODUCTION 

cr 

Typical quantum-mechanical treatments of molecular systems are based on the Born-Oppenheimcr (BO) approxi- 
mation that separates the motions of nuclei and electrons One first solves an eigenvalue equation for the electrons in 
\Q , the field of nuclei clamped at some predetermined points in space and thus builds the potential-energy surface (PES). 

Then one solves an equation for the nuclei moving on that PES and obtains the allowed energies of the molecule. In 
£SJ ■ many cases the calculation is restricted to the neighborhood of the minimum of the PES in order to determine the 
• molecular geometry 2J. Accurate PES's are useful for chemical-kinetic studies^. 

There has recently been great interest on the calculation of molecular properties by means of nonadiabatic ap- 
proaches; i.e. without resorting to the BO approximation. Since the eigenfunctions of the Hamiltonian operator of 
the molecule H are not square integrable with respect to the 3N coordinates that describe the N particles (electrons 
plus nuclei) in the laboratory-fixed coordinate axes one has to remove the unbounded motion of the molecular center 
of mass and the three corresponding coordinates. Thus one is left with the Schrodinger equation for the resulting 
molecular Hamiltonian operator in the molecule-fixed coordinate axes Hm that depends on 3A — 3 coordinates and 
momenta. For brevity, from now on we call those coordinates absolute and relative (or internal), respectively. By 
means of the variational method one obtains an approximate eigenfunction of Hm that should be also eigenfunction of 
the operators that commute with Hm, such as, for example, spin, angular momentum operator in relative coordinates, 



etcflfl. 
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There are several strategies for solving the Schrodinger equation without the BO approximation. One of them 
is based on explicitly correlated Gaussian functions of the relative coordinates of all the electrons and nuclei in 
the molecule (for a comprehensive review see Bubin et alQ). The trial function constructed from those Gaussians 
satisfies the permutational symmetry of identical particles and contains several parameters that are to be determined 
according to the variational method. If the Gaussians are located at one nuclei it is not difficult to choose the 
variational wavefunction to be eigenfunction of the angular-momentum operator. In some cases it is convenient to 
place the centers of the Gaussians on different space points and it is more difficult to force the variational function 
to be eigenfunction of the angular-momentum operator. In these calculations the authors explicitly expressed the 
variational function and the Hamiltonian operator in relative coordinates [6J (and references therein). 

Another widely spread strategy is based on uncorrelated functions and the SCF approach. In this case the variational 
function is written as a product of one-particle functions with the appropriate permutational symmetry. Since there 
are nuclear orbitals in addition to electronic ones this approach has been termed dynamic extended molecular orbital 
(DEMO) methodQ or nuclear orbital plus molecular orbital (NOMO) theory [g- 11 1. Since such orbitals are expressed 
in terms of the 3iV absolute cartesian coordinates one has to be careful to avoid the spurious contribution of the 
kinetic energy of the center of mass 12] . In order to avoid this problem Nakai et al^-ll proposed to subtract the 
kinetic energy of the molecular center of mass from the Hamiltonian operator during the variational optimization of 
the trial wavefunction and developed the translation- free NOMO (TF-NOMO). 

Although the NOMO approach is a nonadiabatic method its implementation is reminiscent of the BO approximation 
in that the orbitals are located at the "experimental geometries" Q . Note that such a concept is rather alien to the 
nonadiabatic quantum-mechanical calculation just outlined (compare it with the more rigorous approach described 
by Bubin et al|6|). The properly symmetrized product of orbitals located at differentpoints and expressed in terms 
of absolute coordinates is not an eigenfunction of the angular- momentum operator |8|- 
rotational states are expected to contribute to the optimized variational function. 

In order to remove the contribution of rotational states with angular-momentum quantum number J > Nakai 
et al[s 11 proposed a rotation- free NOMO that consists of subtracting also the rotational kinetic energy from the 



llj . For that reason several 



total Hamiltonian operator. The resulting approach is named translation- and rotation- free NOMO (TRF-NOMO). 
However, the removal of the spurious rotational kinetic energy in this way is not exact as in the case of the translation 
kinetic energy as argued by SuttclifFejl3|. 

The purpose of this paper is the analysis of the performance of the TF-NOMO and the effect of using absolute 
cartesian coordinates in the trial wavefunction. Since the treatment of realistic examples may be rather cumbersome 
we apply the method to a simple exactly solvable model. 

In order to make this paper sufficiently self-contained, and to facilitate the discussion throughout, in Sec. |H] we 
outline the separation of the kinetic energy of the center of mass and the construction of the molecular Hamiltonian 
in relative coordinates. In Sec. lIIII we apply the TF-NOMO to an exactly solvable model in order to test the accuracy 
of this approach and the effect of using the absolute coordinates instead of the internal ones. Finally, in Sec. IIVI we 
discuss the results and draw conclusions. 
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II. MOLECULAR HAMILTONIAN 

In this section we outline some general properties of the nonrelativistic Hamiltonian operator for a system of N 
charged point particles with only Coulomb interactions. The results are well known and have been discussed by several 
authors in different contexts (see, for example, the review by Fernandez and Echave|5( and the references therein). 
The nonrelativistic Hamiltonian operator for a molecule can be written as 

H = f + V, 



^ 2m, ' 
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where m, is the mass of particle i, Qi — — e or Qi — Z^e are the charges of either an electron or nucleus, respectively, 
and Tij = \vi — Vj\ is the distance between particles i and j located at the points r,; and r^, respectively, from the 
origin of the laboratory- fixed coordinate system. In the coordinate representation p; = —iKVi. 

Since the the uniform translation of all the particles U(a)riU(ay = Vi + a leaves the Coulomb potential invariant 
tl(a)Vtl(ay = V, then the eigenfunctions of the translation-invariant Hamiltonian operator (fTJ) are not square 
intcgrable. For that reason we separate the motion of the center of mass and define translation-invariant internal 
coordinates by means of a linear transformation 

i 

It is our purpose to keep the transformation ([2]) as general as possible so that it applies to a wide variety of nonadiabatic 
approaches. We arbitrarily choose r[ to be the coordinate of the center of mass 

t li = ^,M = Y j m i (3) 

i 

and r'-, j > 1 the translational-invariant coordinates 

£>i< = a, j>i (4) 

i 

Note that if the coefficients of the linear transformation @ satisfy equations (J3]) and ^ then 

L>(a)r;-C/(a)t =r' ] + a8 jx (5) 

The choice of the coefficients of the transformation © for the translational-invariant variables , j > 1 is arbitrary 
as long as they satisfy Eq. ((4]) (for a more detailed discussion see Q). 

As a result of the change of variables, the total Hamiltonian operator reads 

4-2 

" = -2M V ' + "" 

*- = 4ssfe^) ^ ^Y.t 9 ^ (•) 

j>l k>l \ i / u i=l j=i+l 3 
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where Hm is the internal or molecular Hamiltonian operator. The explicit form of the interparticle distances r%j in 
terms of the new coordinates r' k may be rather cumbersome in the general case but there are particular choices that 
are suitable for the calculation of the integrals necessary for the application of the variational method^, 3]. The 
treatment of the simple model in Sec. IIIII shows one of those particular transformations. 

For future reference it is convenient to define the center of mass and relative kinetic-energy operators 

*CM = (7) 
j>l k>l \ i / 

respectively, so that T = T CM + T rel , H M — T rel + V, and H = T CM + H M . 

The inverse transformation t _1 exists and gives us the old coordinates in terms of the new ones: 



According to equations ([5]) and ((9J we have U(a)TiU(a)^ = (t a + rv from which we conclude that 

(t-% =1, i=l,2,...,JV (10) 

In order to understand the meaning of this result note that the momentum conjugate to v\ is given by the transfor- 
mation 

3 

so that the linear momentum of the center of mass 

p'i = E^' ( 12 ) 

3 

is precisely the total linear momentum of the molecule. We also appreciate that 

TCM = ^ = ^]T$>-P; (13) 

»=1 3=1 

The eigenfunctions of the total Hamiltonian operator ([T]) are of the form 

*(r 1 ,...,r N ) = e ik <TP(r' 2 ,...,r' N ) (14) 

with the appropriate permutational symmetry for the electrons and nuclei. We are of course interested in the wave- 
function for the internal degrees of freedom "0( r 2i ■ • ■ ; v 'n) that provides the relevant molecular properties. For this 
reason we should use a trial function of the corresponding coordinates tp(r' 2 , ■ ■ ■ , r^,) and apply the variational method 



with the relative or molecular Hamiltonian operator Hm [ 

Nakai et al[8-ll| proposed an alternative route based on a trial function of the absolute cartesian coordinates 
<ys(ri, . . . , rjv) and applied the variational method to 

w= (<p\H-Tcm\<p) 
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More precisely, they resorted to the Hartree-Fock method with electronic and nuclear orbitals. In this equation one 
integrates with respect to dr\dv 2 ■ ■ . dr^v (including spin if necessary) and the trial function <£>(ri, . . . , r/v) is square 
intcgrable with respect to all those 37V electronic and nuclear variables. This approach is called TF-NOMO method and 
is an improvement over the translation-contaminated NOMO TC-NOMO method (like, for example, the DEMoQ) 
that is based on the variational method for H. Note that the domain of the NOMO trial function y>(ri, . . . , r/v) is 
R 3N and that for the molecular wavefunction ip(r' 2 , . . . , r' N ) is R 3N ~ 3 . Thus, from a quantum-mechanical point of 
view they belong to different state spaces. 

If we rewrite the arguments of the trial function in terms of internal coordinates </?( r ij . . . , rjv) = <^(ri, r 2 , ■ ■ . , r' N ) 
then we appreciate that the probability distribution of the internal variables 

p(r' 2 ,...,r' N ) = J ^(r 1 ,r' 2 ,...,r' N ) 2 dr 1 (16) 

may exhibit some degree of correlation even though the trial function y>(fi, ■ ■ ■ , rjv) is merely a product of nuclear 
and electronic orbitals with the appropriate permutational symmetry [8- 11 1. The analysis of the effect of this particle 
correlation in a realistic molecular system appears to be rather complicated; for this reason in what follows we resort 
to a quite simple example. 



III. SIMPLE MODEL 



In order to have a clearer understanding of the TF-NOMO method [s- 11 1 we apply it to a simple exactly-solvable 
model. We are only interested in the removal of the translational contamination because getting rid of the rotational 
one does not appear to be so simple [3]. Therefore, a one-dimensional model with at least three particles and a 
translation- invariant potential will suffice. Our model consists of three particles of masses Mi, M 2 and M3 that move 
in one dimension and interact through forces that follow Hooke's law: 



H 



h 2 d 2 h 2 d 2 h 2 d 2 



2Mi dX 2 2M 2 dX 2 2M 3 dX 2 



1 

+ 2 



where are the force constants. 



k 12 (x x - x 2 y + k 13 {x x - x 3 y + k 23 (x 2 - x 3 y 



(17) 
(18) 



In order to reduce the number of parameters to a minimum we first assume that the three particles arc identical, 
so that Mi = M 2 = M 3 = M and K\ 2 = K X3 = K 23 — K. We define dimensionless coordinates x, = Xi/L, where 



L = %/{rnui) and ui = y K/M and obtain the dimensionless Hamiltonian 

1 (&_ d 2 . d 2 

~2 \dxj 
(xi 



H 

TlLO 

1 

+ 2 



x 2 ) 2 + (xi - x 3 ) 2 + (x 2 - x 3 ) 2 



(19) 



and the dimensionless energy e = E/(hu>). 

We separate the motion of the center of mass by means of the transformation 



'7i 

'72 
33 



1 / 

- (Xi +x 2 +x 3 ) 



x 2 
x 3 



X\ 



(20) 
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where qi is the coordinate of the center of mass and q 2 and q 3 are the coordinates of particles 2 and 3 with respect to 
the coordinate origin located arbitrarily at particle 1. These q/s are the r^-'s of Sec. HJ and the transformation (|20|) 
satisfies the equations discussed there. We thus obtain 



H A = T, 



CM 



H r 



(21) 



that is the sum of the dimensionless kinetic energy of the center of mass Tcm and the dimensionless Hamiltonian 
operator for the relative motion H r (Hm in the general discussion of Sec. [TTJ) 



1 8 2 

Tcm = --^^n 

6 8q{ 



H r 



if 



dql 



dq 2 3 



The eigenfunctions are of the form 



dq 2 dq: : 



*(9i,«2,«3) = e lkqi ^{ qil q 2 ) 



+ ll + ?3 - 9293 



(22) 



(23) 



where ip(qiiQ2) is an eigenfunction of H r and — oo < k < oo. Note that ^(qi, q 2l 93) is not square integrable with 
respect to dx\dx 2 dx 3 as expected from the fact that the motion of the center of mass is unbounded. The total 
dimensionless energy is ct = ( + fc 2 /6, where e is the dimensionless energy of the relative motion (an eigenvalue of 

H r ). 

If we try an uncorrected gaussian function ip(q 2 ,q 3 ) = [(4a 2 — 6 2 ) /n 2 ] ^ 4 exp [—a {q 2 + q 2 ) — bq 2 q 3 ] , where 4a 2 — 
b 2 > 0, then we obtain the exact ground-state eigenfunction of H r 

V3 



^00(02,03) = ^exp 



(<72 - °20 3 + gf) 



(24) 



with dimensionless energy e o = Vs. 

As uncorrelated variational function in the absolute coordinates we try 

3/4 



f \ f 2a \ 



exp 



4)} 



(25) 



that is square integrable with respect to dx\dx 2 dx 3 . Note that there is a redundant coordinate because we need just 
two variables to describe the bound states of this model as shown in Eq. (|24[) . This trial function is our simple version 
of a NOMO one. The optimal value of the variational parameter a is determined by the minimum of 

(26) 



W(a) = (H~Tcm 



as proposed by Nakai et al(8-lll|. where 

I ( ° 2 
6 \dx 2 



T. 



CM 



d 2 



d 2 



d 2 



dx\ 



dx 2 



dx\dx 2 dxidx 3 



8x 2 8x 3 



(27) 



is the kinetic-energy operator for the center of mass in the absolute coordinates. Note that present Tcm is what 
Nakai et alj^j call Tt and approximate by T n> x in their calculations, and present approach is a simplified version of 
the TF-NOMO. 

At first sight it is surprising that the variational function with the optimal value a = V3/2 of the adjustable 
parameter yields the exact energy W(V3/2) = vo. However, it is not the only striking fact because the expectation 
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values of any function of Xi — X\ and X3 — x\ (like, for example, (xn — xi) 2 , (X2 — £i) 4 , (^2 — 2^1X^3 — %i), etc) are 
exact too. In spite of this remarkable agreement the exact and approximate wavefunctions are not the same as follows 
from the fact that TcMipoo — and I^p\TcmV^ = V / 3/4. We will explain these curious results later on; for the 
time being note that the results of Nakai et al[8j for Ef f c - E? a f = ($ | T„,t |$o) ~ (*&o I Tt |$o) yield the spurious 
contribution of the kinetic energy of the center of mass to the molecular energy. They state that this problem is 
due to the fact that the Gaussian functions are unsuitable for describing the translational energy but it is clear that 
a set of Gaussian functions in relative coordinates will not exhibit such undesirable behavior. In other words, the 
translational contamination is a consequence of adopting the absolute coordinates and not a result of the choice of 
Gaussian functions. The same argument applies to the rotational contamination; Bubin et al[6[ show how to obtain 
Gaussian states with zero angular-momentum quantum number (J = 0). 

We may suspect that the unexpected success of the variational approach is partly due to the symmetry of the 
problem (three identical particles). In order to break it with the slightest modification of our model we choose 
K\2 = K13 = K ^ if 23 and define A = K23/K, so that 

_1 ( JP_ d 2 \ 



2 \dx 2 dx\ dx 2 J 
1 

"2 



{x\ - x 2 ) 2 + (xi - x 3 ) 2 + A (x 2 - x 3 ) 2 (28) 



and 



fl ■ = -|rlr^ + 5 [ ' l+ ' l+Afe -' 3,^1 (29) 

Since the masses remain the same the transformation from absolute to relative coordinates and the form of Tcm are 
still given by equations (l2Tfl) and (|2T|) . respectively. 

The exact ground-state wavefunction and energy are given by 

-0oo(<72,g3) = -j= exp [-a (gf + qf) + bq 2 q 3 ] 



V& (V6A + 3 + 3A + 2) V3A + 2- V6A + 3 

(30) 



12(3A + 1) 
+ 2 - V6A + 3 



and 



\/2(2v / 2ATT+V / 3(A + l)) V3A + 2- V6A + 3 
600 2(3A + 1) (31) 

respectively. 

As in the preceding example we consider an uncorrelated NOMO-like trial function of the absolute coordinates 

tp(x 1 ,X2,x 3 ) = exp \-ax\- b{x\ +x 3 )] (32) 

where the optimal values of a and b minimize W(a, b) = (Hd — Tcm\ In this case we have the following TF-NOMO 
parameters and energy 
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y/3 y/6(X + l) 

2 ' 4 



W(a,b) 



V6(A + 1) 



V3 
3 



(33) 



In order to measure the effect of keeping the kinetic energy of the center of mass we also choose the values of the 
variational parameters from the minimum of W(a, b) = (^Hd\ which yields 

a = ^ 



b = 
W(a,b) = VI 



VA+T 



V2 
2 



(34) 



From now on we refer to equations ([33]) and (|34f as TF-NOMO and TC-NOMO in order to make a connection with 
the approach of Nakai et al 8| . 

Fig. [1] shows the exact and approximate dimensionless ground-state energy calculated in the two ways just outlined. 
As expected, TF-NOMO yields considerably more accurate results because TC-NOMO is strongly contaminated with 



12J with respect to the DEMO 



the kinetic energy of the center of mass. This point was discussed earlier by Fernandez 

n 

method 7]. Note that TF-NOMO is remarkably accurate for all A and yields the exact result for A = 1 as discussed 
above. Besides, the approximate energy is smaller than the exact one for A > 1 reflecting the fact that the variational 
principle does not apply in this case because the exact and approximate functions belong to different spaces. 

We may try and improve the TC-NOMO results by simply subtracting the kinetic energy of the center of mass, 
thus producing a sort of corrected TC-NOMO or CTC-NOMO: 



W 



CTCNOMO 



TCNOMO\ 



Hd — Tqm <p 



TCNOMO 



) 



(35) 



Figure Q] shows that the energy calculated in this way agrees quite well with the exact and TF-NOMO ones. This 
result shows that most of the error in the energy calculated by means of the TC-NOMO comes from the spurious 
kinetic energy and just a relatively small contribution comes from the inadequate optimization of the variational 
wavefunction with respect to Hd- 

It is not difficult to explain why the uncorrclatcd trial function in absolute coordinates yields such good results 
(even exact ones for A = 1). It we substitute X2 — <Z2 + x\ and £3 = q% + x\ into equation (|32p and integrate the 
square of ip(xx, q2, 93) = <f(xi, q^ + x%, 93 + Xi), with respect to x\, we obtain 



^(92,93) 



(p(xi,q 2 ,q3) 2 dxi 



2^/ab 



: exp 



2c (a + b) . 2 
' a + 2b [Q2 



4& 2 g 2 g3 
a + 2b 



(36) 



TiVa + 2b 

Note that the use of the absolute coordinates in the trial wavefunction introduces some sort of correlation between the 
translation-invariant coordinates qi and 03 when we integrate with respect to the redundant variable. The resulting 
correlation is reasonable because it is determined by the variational method, and, in particular, when A — 1 it yields 
(fortuitously) the exact probability distribution for the relative coordinates 



p{q2,qz) = - exp 

7T 



^00(02, q?,Y 



(37) 
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It is now clear why we obtained the exact energy and expectation values before for this particular case. In fact, we 

expect to obtain the exact expectation values of any observable in relative coordinates. We do not obtain the exact 

expectation value of Tcm because this operator contains a derivative with respect to the redundant absolute variable 

that does not appear in the exact square-integrable wavefunction. 

In order to compare the approximate and exact wavefunctions we write both ^(92,93) and "000 (92, 93) 2 as 
1/2 1 

(4a 2 — f3 2 ) /(2"7r) exp [—a (q| + g|) + ^9293] (2a > \f3\). Figures [5] and [3] show that the exponential coefficients 
a and ft, respectively, given by the TF-NOMO agree remarkably well with the exact ones, whereas the TC-NOMO 
coefficients are considerably less accurate. Note that the effect of keeping the kinetic energy of the center of mass 
not only affects the energy (which is expected) but also the form of the variational wavefunction. The agreement 
between the exact and TF-NOMO exponential parameters a and /3 explains why the approximate wavefunction yields 
so accurate expectation values of operators in the relative coordinates qi and q^. In particular, the TF-NOMO results 
are exact for A = 1, but, of course, such a result is not to be expected in a realistic case. 
If we minimize ((p\ H r \(p) with an uncorrelated trial function of the relative coordinates 

<^(<?2, 93) = y^ ex P [-a(<?2 + iD] ( 38 ) 



we obtain the optimal variational parameter a — y/2X + 2/4 and the resulting approximate energy W(a) = \/2X + 2 is 
considerably less accurate than that given by Eq. (|33|). This result suggests that it is preferable to apply the NOMO 
method with orbitals that depend on the absolute coordinates as long as we remove the kinetic energy of the center 
of mass in the optimization process. That is to say, we minimize the expectation value of H r = H — Tcm expressed, 
for simplicity, in the same set of absolute coordinates chosen for the NOMO variational wavefunction. Under such 
conditions the TF-NOMO-SCF method appears to take into account part of the correlation energy that in the present 
simple model is given by W(a,b) — W(a). This energy difference is quite similar for TF-NOMO and CTC-NOMO 
according to Fig. [TJ 



IV. CONCLUSIONS 



We have tried to elucidate the effect of using absolute coordinates in the NOMO-SCF variational method. Since 
such an analysis for an actual molecule, even as simple as H2, is rather complicated we chose a simple model of 
three particles with harmonic interactions in one dimension. Although rather oversimplified, this model enables us to 
take into account the main ingredients of the NOMO-SCF method. We have a translation-invariant potential-energy 
function and, consequently, we have to remove the unbounded motion of the center of mass. In this case a correlated 
Gaussian function of the relative coordinates yields the exact result that is most convenient to test the approximate 
ones. 

The NOMO-SCF wavefunction is simply a product of Gaussian functions (orbitals) for each of the particles. The 
integration of the square this uncorrelated function with respect to the redundant absolute coordinate (three in a 
realistic case as shown in Eq. (TTBIO gives rise to some kind of particle correlation. If the adjustable parameters in this 
trial function are optimized variationally with the Hamiltonian H r = H — Tcm then the correlation just mentioned 
appears to improve the calculation of the energy and expectation values considerably. Perhaps, one should not expect 
such a remarkable success for an actual molecule, but, however, there is no doubt that even in that case the NOMO- 
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SCF will take into account part of the correlation energy in spite of being based on uncorrelated Gaussian functions. 
In order to verify this conjecture that is expressed in Eq. (|16j) it is only necessary to calculate the energy with an 
uncorrelated NOMO trial function and H r both in terms of relative coordinates (like present Eq. (|38|l). Note that 
this interesting feature of the TF-NOMO has apparently passed unnoticed in the applications of the method^- 11 1. 
This fact reinforces our claim on the utility of simple models for the study of rather complicated problems. 

The simple model also shows that if we simply subtract the expectation value of the kinetic energy of the center of 
mass then the resulting energy is quite accurate (what we have called CTC-NOMO). However, in order to improve 
the calculation of the expectation values of other observables it is convenient to apply the SCF procedure with the 
relative Hamiltonian operator H r = H — Tcm that leads to what is commonly known as TF-NOMO [g- 11 



Comparison of the exact and approximate energies for the simple model discussed in the preceding section also 
enabled us to verify that the TF-NOMO may not always yield upper bounds because the variational trial function of 
the absolute cartesian coordinates and the exact eigenfunctions belong to different state spaces. For this reason the 
variational principle does not apply to the TF-NOMO even though this approach is based on the HF-SCF method. 
This unwanted feature of the TF-NOMO does not appear to be a serious drawback as shown by the accurate results 
obtained in Sec. 
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FIG. 1: Exact (solid line, black), TF-NOMO (dashed line, red), TC-NOMO (points, blue) and CTC-NOMO (circles, green) 
ground-state dimensionless energy 
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FIG. 2: Exact (solid line, black), TF-NOMO (dashed line, red) and TC-NOMO (points, blue) exponential coefficient a 




FIG. 3: Exact (solid line, black), TF-NOMO (dashed line, red) and TC-NOMO (points, blue) exponential coefficient /3 



